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Utilizing a number of results of Dittmann, we investigate the nature of the Yang-Mills field over 
the eight-dimensional convex set, endowed with the Bures metric, of three-level quantum systems. 
Parallelling the decompositions of eight-dimensional Euclidean fields by Corrigan, Devchand, Fairlie 
and Nuyts, as well as Figueroa-O'Farrill and others, we investigate the properties of self-dual 
and anti-self-dual (fi-) four-forms corresponding specifically to our Bures/non-Euclidean context. 
For any of a number of (nondegenerate) 3x3 density matrices, we are able to solve the pair of 
eigenequations, *F = ±(^)(n± A F), where * is the Hodge operator with respect to the Bures 
metric and F a two- form. The resultant sets of (traceless) twenty-eight A+'s coincide with the sets 
of twenty-eight A_'s, the sets consisting of four singlets and three octets. The four-forms are 
found to exhibit quite simple behaviors, though we are not able to derive them in full generality.. 

PACS Numbers 11.15.-q, 03.65.-w, 02.40. Ky 



The Bures metric, defined on the nondegenerate density matrices, has been the object of considerable study [Q-^. 
It is the minimal member of the nondenumerable family of monotone metrics . Other members of this family of 
particular note are the "Bogoliubov-Kubo-Mori" (BKM) metric |^,||, the maximal monotone metric, as well as the 
"Morozova- Chentsov" Q and "quasi-Bures" pT] | ones, the last yielding the minimax/maximin asymptotic redundancy 
in universal quantum coding [ p^ . Interestingly, these (operator) monotone metrics correspond in a direct fashion to 
I certain "measures of central tendency" , with, for example, the Bures metric corresponding to the arithmetic mean, 
{x + y)/2, of numbers x and y |^;|^. All these monotone metrics constitute various extensions to the quantum domain 
of the (unique) Fisher information metric in the classical realm where the objects of study are probability distributions 
(rather than density matrices) (cf. |ljl). However, the only one of these monotone metrics that can be extended 
to the boundary of pure states yielding the standard Fubini-Study metric on this boundary is the Bures (minimal 
monotone) one sec. IV] (cf. ||l5[] ). 
I Dittmann has shown that "the connection form (gauge field) related to the generalization of the Berry phase to the 
mixed states proposed by Uhlmann satisfies the source-free Yang-Mills equation *D * Dlu, where the Hodge operator 
is taken with respect to the Bures metric on the space of finite-dimensional density matrices" [|6| (cf. jl^]). "These 
findings may be seen as extensions to mixed states of numerous examples relating the original Berry phase to Dirac 
monopoles, and the Wilczek and Zee phase" [|8|. (Let us also indicate here that the BKM metric has been shown 
to be the unique monotone Riemannian metric for which two "natural flat" connections — the "exponential" and 
[ "mixture" ones — are dual |^.) 

In sec. 6], Figueroa-O'Farrill studied the eigenspaces of various endomorphisms (T), defined by four-forms (T) 
in E^, of the space of two-forms. ("One can use T to define a generalised self-duality in eight dimensions by demanding 
that a two- form belong to a definite G-submodule of bos. This generalises self-duality in four dimensions, where we 
can take T = *1, and T = * itself. The eigenspaces of T are the subspaces of self-dual and anti-self-dual two- forms" 
pof.) We seek to extend this type of analysis to our different (non-Euclidean/Bures), but still eight-dimensional 
setting. 

To proceed, we exploit our recent work pO| , pl| in determining the elements of the Bures metric using a certain 
Euler angle parameterization of the 3x3 density matrices (Earlier, Dittmann had noted "that in affine 

coordinates (e. g. using the Pauli matrices for n = 2) the [Bures] metric becomes very complicated for n > 2 and 
no good parameterization seems to be available for general n".) This allows us, among other things, to employ as 
our parameter space four (noncontiguous) eight-dimensional hyperrcctangles (each of the four having three sides of 
length TT, three sides of length ^, one j in length and one cos~^ "Tj)' ^^^her than less analytically convenient ones, 
such as that discussed by Bloore Fig. 3] , involving spheroids and "tetrapaks" . (The four hyperrcctangles are not 
adjacent due to the fact that two of the six Euler angles used in the parameterization have disconnected ranges. This 
situation, that is a multiplicity of hyperrcctangles — which leads to different normalization factors only — has only 
come to our attention relatively recently, and serves as an erratum to earlier analyses (24|j2^ , which had relied upon 
work of Marinov |2^, without taking into account a correction |2^] to ||2^] (cf. [p8|).) 

For the required connection form A, Dittmann has presented the general formula ]^, eq. (9)], 
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A= ^^^—.{W*T -T*W). (1) 

L -\- R 

Here the elements of T lie in the tangent space to the principal [/(7i) -bundle, the manifold of invertible normalized 
(TrT/F*VF = 1) Hilbert-Schmidt operators, while L is the operator (depending on W) of left multiplication by the 
density matrix p = WW* and R the corresponding operator of right multiplication. Also, L and R are the counterpart 
operators for p = W*W. Since the Euler-angle parameterization of the 3x3 density matrices presented in is of 
the ("Schur-Schatten") form, 

p = UDU\ (2) 

U being unitary, U* its conjugate transpose, and D the diagonal matrix composed of the three eigenvalues of p, one 
can immediately express W as UD^/'^U* , with W = W*. Let us note here that Sjoqvist et al have addressed the issue 
of mixed state holonomy — first raised by Uhlmann (l) as a "purely mathematical" problem — developing "a new 
formalism of geometric phase for mixed states in the experimental context of quantum interferometry" |]30| . In their 
notation, 

n = iTrlpoW^'dW] (3) 
is regarded as a gauge potential on the space of density operators, where 

the unitary operator U{t) being used to fix the parallel transport conditions. It would be of interest to examine 
the relation of their expressions to that of Dittmann (|l|), which we implement here (cf. [^). (In the context of 
pure states, Sanders, de Guise, Bartlett, and Zhang have presented a scheme for producing and measuring an Abelian 
geometric phase shift in a f/iree-level quantum system where states are invariant under a non- Abelian group . They 
geodesically evolve C/(2)-invariant states in a four-dimensional SU{3)/U{2) space, using a three-channel interferometer 
(cf. II).) 

In the non-trivial task (involving computations with modular operators [Q) of implementing formula for the 
three-level quantum systems, we relied upon the implicit relation (equivalence) between two formulas for the Bures 
metric for n-level quantum systems p^ , eqs. (2) and (16)], 



and 



n 

(The somewhat involved formula that we used for the coefficients aij , functions of elementary invariants, is stated in 
Proposition 3 of |^^.) In other words, we used the formula for the connection 
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A = -J2a,,p^-'Sp^~\ (7) 

where S = W*T - T*W . 

For n = 3 we have determined, for particular points in the eight-dimensional manifold, four-forms $1+ for which 
= 17+, where the Hodge star is with respect to the Bures metric as well as four-forms il- for which 

*r2_ = ri_. We formulated this problem as a set of seventy linear simultaneous homogeneous equations in seventy 
unknowns (that is, unknown constants). Selecting various points in our eight-dimensional space (and making use of 
exact arithmetic) , MATHEMATICA has consistently succeeded in expressing thirty- five of the unknowns in terms of 
the other thirty-five variables. However, we did not find a a unique such solution valid for arbitrary points in the 
eight-dimensional space, with MATHEMATICA apparently indicating that such a general solution over the manifold 
with constant coefficients for il is not possible. We have obtained precisely the same type of findings in our search for 
anii-self-dual four-forms. (The Hodge * operation squares to 1 on four-forms in eight dimensions. It has eigenvalues 
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+1 and -1 with equal multiplicity. So one can always find a basis where * is diagonal with thirty-five -|-l's and thirty- 
five -I's.) In all these calculations we have used the duality relation, e g. (1.5.15)] (making use of our explicit — 
though, in some cases somewhat cumbersome — formulas for the g'^ 's [ |20| , ^ ) 



where we took n = p = 4 

For the specific point qi, we set the eight variables, described in |2^, tobea = — ^,t = ^,a = ^,6 = — 



^,9 = —^,01 = —^,02 = — f — the first six of these being Euler angles that parameterize the unitary matrix U, 
while the last two parameterize the diagonal matrix of eigenvalues D in the decomposition (^. In our eight-dimensional 
manifold of states, the point qi gives us the specific 3x3 density matrix, 

/ 307 3(59i+25\/3) (-l)^/^^(51+29ix/3) 

' 1024 2048 1024 

3(-59i+25V3) 417 3(-l)^^^''(2li+25\/3) (9) 



Pi 



2048 1024 1024 

V5T2(-l)"/''(9 + 20iV3) ^(-1)3/4(24 + iV3) ^ 



This matrix is, of course, Hermitian, nonnegative definite — possessing eigenvalues — and has trace 1. (It 

corresponds to a strictly impure or mixed state. Also, we should point out that the entries of the 8x8 Bures metric 
tensor turn out, in fact, to be independent of the variables a and a, as established in ]20|.) For the solutions for A+ 
at qi of the eigenequation, (F being a two- form), 

*F=^in+AF), (10) 

(having substituted unity for each of the free thirty-five variables in the thirty-five-dimensional solution to the set of 
seventy linear homogeneous equations for 17, based upon (^, which enforce self- duality), we obtained the twenty-eight 
eigenvalues (cf. eq. (2.1)]) 

X+{qi) 6.15149, -6.06045, ±5.11128(eightfold), -4.16211, 4.07107, ±.994689(eightfold), ±.0455182(eightfold). (11) 

Of course, the three octets are each split into two quartets of opposite sign. (The adjoint representation of 5[/(3) is 
eight-dimensional, with its elements generically lying in 50(8) sec. 7].) The four singlets or isolated values are 
the roots of the quartic (biquadratic) equation, 

17848517231861271296 + 718875A+(-72767012864 - 2131611323040A+ -I- 39304490625A3_) = 0. (12) 

The six possible values assumed by the three octets are the roots of the sextic equation (a cubic equation in x^), 

82734971267961585664- 18225A^(2195802859754043904+ 291144375A^(-7894856752 -I- 291144375A^)) = 0. (13) 

(Consistent with tracelessness, the twenty-eight eigenvalues sum to zero — as well as, therefore, the four isolated 
eigenvalues themselves.) We have (as one particular example) "simplified" the exact symbolic expression for the root 
of ( |l^ ) corresponding to the eigenvalue 5.11128, so that it is exphcitly real- valued, obtaining thereby that 

(5.11128)2 w ( 16( 493428547+ (14) 

^ ' 873433125V V ^ ' 



2V217739666231788507 cos 



19986057V257834787813597115559383045701069731 



101094855629270248323646732 



(The value of the cosine term here is close to unity, that is .999444. Any further simplification of (Q) does not seem 
possible.) 

We should also note that the cardinality of the set of eigenvalues, that is twenty-eight, corresponds to the number of 



entries of the (antisymmetric) two-form F, that is Fij = 1, 



), for which i < j (having observed that Fj,, 



-F- ■ 

and Fa = 0). Proceeding in the exact same manner, but using the anti-self-dual form MATHEMATICA provided, 
and solving 
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*F = -^in-AF), (15) 



we obtained precisely the same set of twenty-eight eigenvalues (|l 1\) , despite that fact that Q+ ^ ±i7_. 

We observe here that the well-known octonionic instanton (CDFN) equations are obtained by choosing a 
particular constant (hence closed) four- form in R*. This form decomposes the bundle of two- forms into two sub- 
bundles of dimensions seven and twenty-one. Each of these sub-bundles corresponds to a different nonzero eigenvalue 
(-3 and 1, respectively) of the transformation (which is symmetric, hence diagonalizable, as well as traceless) 

F^~^^{nAF), (16) 

so restricting a Lie-algebra valued two-form to lie in either of these two sub-bundles guarantees that it will satisfy the 
Yang-Mills equations of motion. 

We have (denoting by Qjki the four form dxi A dxj A dxk A dxi, where for our eight variables we take xi = a,X2 = 
T,X3 = a,X4 = (3,x^ = 6, = OjXi — 9i,xs — 02, and using lexicographic ordering of the indices), corresponding to 
the point qi, 

^ , , 378375 . 14037 . 601 59079 . . . 

— ig54ox6 ^^^'^ ~ 127232 " ' ^^Ci578 284~ ''^ ~^ '^^ " ' ^^^"^^^ ^ ' 

and 

6975 14187 647 58745 

i'-(9i) — 23296 ~ 127232 ^^^'^ ' " — 7]~''i578 + 284 ~^ '^^ " ' ~^ Csers- (18) 

So we see that ^2+(gi) ^ ±^-{qi) (nevertheless, clearly sharing certain features). When we used these same two 
four-forms, but at another point (a = r = y^, a = = = j,9 = |,6'i = |,02 = f) than that yielding pi, 
the two sets of eigenvalues for the equations and (|l|) were quite similar (e. g. listing the largest ones in absolute 
value, 9.83657 vs. 9.66359, -9.73817 vs. -9.59167,...) but not now strictly equal, and the twenty-eight eigenvalues in 
each set were all distinct from one another. 

We, then, studied a second density matrix (cf. (0)), 



128 32 ^ 128 ^ ' ' \ J 



P2 



32 ' 128 
41 

32 128 128 128 128 



^_151 AL ^_^_^)((„l)5/12 + (_i)3/4) 



(19) 



V (lis + ii)((-i)'/' + ^^+1^^-^^+-^^ i J 



having eigenvalues 5, | and i. It corresponds to a new point q2 for which 

TT 3tT 27r TT TT TT TT TT , , 

a^-,r^-,a^-,b^-,P=-,9=-,e, = -,e2 = -. (20) 

Adopting the same analytical strategy as employed for pi , we obtained for the solutions for A+ at q2 of the equation 
(|l0|), the twenty-eight eigenvalues (cf. (|ll|)), 

X+iq2) = 2.68934, -2.60397, ±2.14857(eightfold), -1.69317, 1.6078, ±.498082(eightfold),.0426838(eightfold). (21) 

So again, fully analogously to (0), we have four singlets and three octets. (Tracelessness is again verifiable. Let us 
also note that for the analysis based on pi, as well as that on p2, the four singlet eigenvalues have the largest, second 
largest, and fourth and fifth largest of the seven possible distinct absolute values.) Computations of the analogues 
based on p2 of the exact results (^2|), ( |l3|) and (14) proved now to be much more intractable, though. The self-dual 
four-form here employed was 

448+ 128a/3 + 27V6 , -3 + 224^^ 

^'+W2j — <,1234 • • • H Q <,1678 + ^,2345 ■ • • + 1,5678: l^^J 

while the anti-self-dual four-form was 

^ , , 448+ 128\/3- 27\/6^ -3 - 224^/6^ ^ ^ 

"-(?2) — ggg C1234 ■ • • H g Cl678 + C2345 . . . + C5678- (23) 
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We found the same set of twenty-eight eigenvalues (21) relying upon ri_((72), using the equation (p^, as we did on 
r2+(g2), using the equation ([Tc|). 

For several additional density matrices (other than pi and p2), we have been able to obtain fully analogous results. In 
all these instances, we let the "second" (lexicographically speaking) set of thirty-five variables (that is, the coefficients 
of C2345 • • ■ Csers) be unity, and solved for the coefficients of the "first" set of thirty- five (C1234 ■ • ■ Ciers)- 

We have also been able to derive univariate generalizations of the (constant) four- forms (22) and (23), allowing 
either the Euler angle r, (3, h or 9 to be free, rather than fixed at particular values. It turns out — as numerical 
evidence firmly convinced us — that the twenty-eight eigenvalues ( pl| ) are completely unchanged as [3 is varied (even 
though the coefficients of the extended four-form themselves do explicitly depend on /?), but do change with t, h or 
9. For example, for the point q2 (pOj) we have obtained the four-form (22). The coefficient there of C1234 is 



C1234 



448-1- 128\/3-F 27^6 
896 



.821249. 



(24) 



If we let the Euler angle r vary (its full range being [0,7r]) rather than being fixed at 37r/4 (as it is, of course, for 52 
itself), we obtain instead of (p3) the more general expression. 



C1234(t) 



(25) 



28 V2(42 - 5\/3) cos 2t + 1 19^6 cos 4t - 4480(7 + 2^3) sin 2r + x/6(2009 - 10 sin 4t) + 84^2 sin 4t 



62720 



Similarly, letting /3, b and 9 vary in turn we derive. 



Ci234(/3) = ;-i-(448-f 128%/3 + 27V6)sin2/3, 
995 



(26) 



Cl234(&) 



1 V3(128 27\/2) sin 2b 

2 ^ 128(7 cos 26) ' 



(27) 



and 



01234(6*) cos6' 



2 sin 6* 9cos6'sin-^6' 
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28^2 



(28) 



Also, 



Cl234(^l) = (29) 

1600 + 256(7 + 8^3) cos 29i -f 64(1 1 14V3) cos A9i + 3\/3(384 -I- 26 sin 9i + 35 sin ?,9i + 25 sin 56*1) 

128(25 -f 28 cos 26(1 + llcos46'i 

(We have not so far — due to increased computational burdens — been able to similarly analyze the scenarios in 
which 02 is free, nor in which pairs of the parameters can simultaneously vary.) In Figures 1 - 5 we plot these five 
simple well-behaved univariate functions (|2^)-(|9|), which generalize the constant coefficient (24) over the ranges of 
their respective variables. 
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C1234 (t) 



. 8 
. 6 
. 4 
0.2 




FIG. 1. Univariate generalization (p5|), incorporating the Euler angle r, of the constant coefficient ( p^ of ^1234 in the self-dual 
four- form fl+{q2) (p^, corresponding to the density matrix p2 ( p9| ) 



C1234 (/3) 




/3 



0.25 0.5 0.75 1 1.25 1.5 



FIG. 2. Univariate generalization (pq), incorporating the Euler angle /3, of the constant coefficient ( |24| ) of 1^1234 in the self-dual 
four-form f2+(g2) (P^), corresponding to the density matrix p2 ([LSJ) 
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C1234 (e) 




FIG. 4. Univariate generalization (|28|), incorporating the Euler angle 9, of the constant coefficient ( |24[ ) of ("1234 in the self-dual 
four- form ^+(52) (p4)' corresponding to the density matrix p2 (|l9|) 



C1234 (01 ) 



1 . 3 



1.2- 



1 . 1 




FIG. 5. Univariate generalization (j29|), incorporating the angle 6^1, of the constant coefficient ( p^ of (^1234 in the self-dual 
four-form Q,+ {q2) (p^, corresponding to the density matrix p2 ( |l9| ) 

The maximum of Fig. || is precisely the value (^), that is .821249. (We already noted that the twenty-eight 
eigenvalues associated with the point (72 are unchanged as /3 varies, in contrast with the other four variables utilized 
here.) The maxima for the other four figures are r = .821472, b = .83522, d = 1.04598, and di = 1.38035, all of which 
exceed .821249. This last value thus falls (as it must) within the ranges of values assumed by the functions above of 
the form ci234(). (Let us also remark that the specific values of the eight parameters employed to derive pi do not all 
strictly fall within the relatively narrow ranges — devised so as to avoid duplication — designated for them in . 
We have so far been unable to fully determine the corresponding eight parameters, within these designated ranges, 
yielding pi. The point corresponding to p2, on the other hand, does lie within the domain of parameters specified in 

As a futher illustration of the apparent simple characteristics of the self-dual four- forms in our analytical context, 
we give the counterparts of (|24|)-(p9|) and Figs. 1-5, but based on the coefficient of Ci678 rather than of C1234, 



ci678(r) = 112^1+ (i + 16V2) cos 2r -I- cos r sin r, (31) 

C1678 (/3) = cos 2/? i (-3 + 224^6) sin 2/3, (32) 
6 
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Also, 



1 2 

C1678 (&) = - 2 + 56 y - CSC 5 sec 6, 



C1678 {9) = - COS 6* + 28 V2 esc 6* sec 6*. 



1 16(13 sin 6ii + sin 36*1) 
£1678(^1) = -7: + 



2 x/3(l + 7cos26ii) 



(33) 
(34) 

(35) 



C1678 (t) 



110 



100 



90 



80 




0.5 1 irs"^ 2 2.5 3 

FIG. 6. Univariate generalization (^), incorporating the Euler angle r, of the constant coefficient ( |30| ) of ^lers in the self-dual 
four-form f2+(g2) (p2|), corresponding to the density matrix p2 (|l9|) 



C1678 (/3) 




/3 

0.25 0.5 0.75 1 1.25 1.5 
FIG. 7. Univariate generalization (^), incorporating the Euler angle of the constant coefficient (|3^) of ("lers in the self-dual 
four- form f2+(Q2) (p^, corresponding to the density matrix p2 ( p^ 
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C1678 (b) 



6000 
5000 
4000 
3000 
2000 
1000 



0.25 0.5 0.75 1 1.25 1.5 
FIG. 8. Univariate generalization (^), incorporating the Euler angle b, of the constant coefficient (^o|) of ^lers in the self-dual 
four- form fl+{q2) (^3), corresponding to the density matrix P2(|l9[) 



Cl678 

5000 
4000 
3000 
2000 
1000 



(9) 



0.25 0.5 0.75 1 1.25 1.5 
FIG. 9. Univariate generalization (p^), incorporating the Euler angle 9, of the constant coefficient ( po| ) of ("lerg in the self-dual 
four-form Q,+ {q2) (Py), corresponding to the density matrix p2 ([LSJ) 



Cl678 

15000 
12500 
10000 
7500 
5000 
2500 



(Si) 



0.2 0.4 0.6 0.8 

FIG. 10. Univariate generalization (^5|), incorporating the angle 9i of the constant coefficient (|3C|) of ^1678 in the self-dual 
four- form Q,+ {q2) (P4), corresponding to the density matrix p2 (u 
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In all this line of work, we hope to further simplify the existing formulas [gOj] for the elements of the Bures metric 
over the 3x3 density matrix and of its inverse, so that the challenging computations of the type reported above might 
be facilitated and further advanced to include the simultaneous consideration of additional variables. 

By way of comparison/contrast of our results, let us note that for the well-known (self-dual) "Cay ley calibration" 




T — C1234 + Cl258 — Cl267 + C1357 + Cl368 " Cl456 + Cl478 (36) 



+C2356 — C2378 + C2457 + C2468 ~ C3458 + C3467 + C5678 , 

which is invariant under a Spin^ subgroup of SOg, the 28 or adjoint representation sos of SOg breaks up as 

28 ^ I ® 21, (37) 

where the 21 corresponds to the adjoint representation spin7 G sog- The endomorphism T of the space of two-forms 
(F) obeys the characteristic polynomial, 

(t -ij(T + sij = 0. (38) 

The eigenvalues are, therefore, 1 and -3, and (consistent with tracelessness) have multiplicities 21 and 7, respectively. 
"Therefore there are two possible extensions of self-duality, and hence two possible extensions of the notion of instanton 
to eight dimensions" (The analogous characteristic polynomials of a number of other possible generalized 

self-dualities — complex, special lagrangian, complex lagrangian, quaternionic, and sub-quaternionic — in eight- 
dimensional Euclidean space are also given in ||l^, sec. 6].) Of course, it would be desirable to develop explanations 
comparable to these of the results we have presented above. 
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